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Abstract. In this talk I'll review the present status of charged particle multiplicity 
measurements from heavy-ion collisions. The characteristic features of multiplicity distributions 
obtained in Au-|-Au collisions will be discussed in terms of collision centrality and energy and 
compared to those of p+p collisions. Multiplicity measurements of d-|-Au collisions at 200 GeV 
nucleon-nucleon center-of-mass energy will also be discussed. The results will be compared to 
various theoretical models and simple scaling properties of the data will be identified. 



1. Introduction 

The multiplicity of charged particles emitted in a heavy-ion collision constitutes an important 
observable, which reflects the properties of the hot and dense system formed in the overlap region 
between the two incoming nuclei. Even without more detailed and differentiated measurements 
of the emitted particles one can obtain important information about the collision from 
measurements of the total multiplicity of charged particles, its distribution in pseudorapidity 
space (angular dependence) and its dependence on collision centrality and energy. Central 
questions concerning the redistribution of the incoming energy into particle production and 
kinetic energy can be addressed on the basis of relatively simple multiplicity measurements. 
It is therefore very appropriate that the topic of particle multiplicity in heavy-ion collisions is 
addressed in this workshop, and in this article I will attempt to give an overview of multiplicity 
measurements at RHIC energies and to compare these data with measurements at lower energies 
as well as simpler collisions between individual protons and leptons. 

Some of the general properties of charged particle multiplicity distributions may be seen 
in Fig. n Here the distribution of charged particles, dNch/dr], where r] = — lntan(0/2), is 
shown as a function of pseudorapidity for central Au + Au collisions P (solid points). The 
solid curve represents a triple-Gaussian fit to the data. One observes several distinct features 
of this distribution: an approximately fiat region near mid-rapidity extends in this case out to 
\t]\ ~ 2 followed by a smooth fall-off region toward larger /smaller values of rj. Since the particles 
emitted within the angular region 45° < 9 < 135° corresponding to —0.88 < rj < 0.88 (grey 
bands in panels a) and b)) are most likely to represent a thermalized region of phase space, 
there is special significance attached to the number of charged particles emitted in this region, 
i.e. the height of the mid-rapidity plateau or dNch{\v\ < l)/dr]. Although this region appears 
as a plateau in the ciA'^c/i/dTy-distribution shown in panel a) it is interesting to observe that the 
dNf-h/dO and the dN^h/dfl distributions both exhibit a minimum at 9 = 90° corresponding to 
r] = 0. In reality the charged particle distributions are very forward peaked. 




Figure 1. Illustration charged particle distribution for 0-6% central 200 GeV Au+Au 
collisions 1 . Panel a: The pseudorapidity density, dN/dr] is shown as a function of rj. Panel b: 
The corresponding angular distribution dN/dO is shown as a function of the angle 9 relative to 
the beam axis. Panel c: same as for panel (b) but here dN/dQ is shown. The shaded regions 
in panels a) and b) indicate the angular region where the transverse momentum pt exceeds the 
longitudinal momentum 



2. Mid-rapidity densities 

Before the start of the RHIC program, several theoretical predictions had been made concerning 
the density of charged particles at mid-rapidity, dN/dy\y=o for ^/sJfN=200 GeV central Au+Au- 
collisions. These predictions are summarized in the left panel of Fig. [3 which has been adapted 
from a compilation by Eskola |2]. The vertical band represents the PHOBOS measurement of 
(iA^/ci?7||^l<i multiplied by a factor of 1.1 2^ to correct for the transformation to rapidity, y. It is 
evident that most of the predictions overestimated the density by up to a factor of two although 
a few predictions agree with the measurement. 

Measurements of dNch/drj may also be used to obtain an estimate of the energy density 
achieved in the collision region. A simple estimate may be based on the naive assumption 
that the total available energy in 200 GeV Au + Au collisions of ~39 TeV could be converted 
into a completely stopped source of high energy density that emits particles isotropically in the 
final state. This leads to an energy density of ^ 130 GeV/fm^, based on a formation time of 
To ~ 1 fm/c. Since this estimate exceeds the critical energy density for the transition into a 
partonic state (QGP) by two orders of magnitude, one may consider it inevitable that this new 
state of matter is formed in such collisions. Studies at lower collision energies performed at the 
AGS and SPS have shown, however, that this complete stopping picture is invalid. Also the 
present measurements are at variance with this picture. The dNch/di] distributions observed are 
substantially wider than the dNch/dr] oc 1/ cosher/ shape corresponding to the fully stopped / 
isotropic source scenario. Alternatively, the energy density may be estimated more reliably from 
the method proposed by Bjorken [7j, which is based on the total energy of particles emitted at 
mid-rapidity. Espressed in terms of the dN^h / drj distribution we find 

_ {mt) dNch ^ . , . 

^0 — , JneutJ Jaci 

where {rrit) is the mean transverse mass, vri?^ is the transverse area of the fireball, tq is 
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Figure 2. Left panel: A compilation of theoretical model predictions (solid squares) of the mid- 
rapidity density of charged particles dNch/dij or dNch/dy made prior to the start of the RHIC 
in July 2000 121. The vertical line is then PHOBOS measurement |3| at 200 GeV multiplied by 
a factor of 1.1 to correct for the dN/dr} to dN/dy transformation. Right panel: The Energy 
dependence of the mid-rapidity |?7| < 1 pseudo-rapidity density. Data obtained at RHIC (solid 
circles) [30]; SPS (open squares) JB], and AGS (open diamonds) 5 show a roughly logarithmic 
increase with collisions energy. y/s^N (solid line) , but are also consistent with the trend of the 
Gluon Saturation Model (dashed curve) at y/SNN > 20 GeV. 

formation time and the factors fneut and fjac correct for the unobserved neutral particles and 
the Jacobian in the pseudorapidity to rapidity transformation. Inserting values of (mt)=0.57 
GeV/c^, /„eMt=l-6, /jac=l-l, all of which are obtained from particle identified spectra measured 
by BRAHMS using dNch/dr]=700 for the 3% most central Au+Au collisions and a formation 
time of To=l fm/c one finds an energy density of eo ~5 GeV/fm'^ , which is still substantially 
above the predicted threshold for the phase transition of eo ~0. 7-1.0 GeV/fm'^ obtained from 
lattice QCD calculations i9j . Much of the available energy is thus carried off by particles 
emitted in forward-backward angles associated with lower energy density, which gives rise to 
pseudorapidity distributions that are substantially broader than expected for an isotropic source. 
The pseudorapidity distribution is therefore intimately connected to the energy density of the 
emitting source and provides an important test-bed for validating or discrediting theoretical 
models attempting to describe the conditions in the early phases of the collision. 

The dependence on collision energy is shown in the right panel of Fig. [21 The mid-rapidity 
charged particle density dNchiv = 0)/dr]/{Npart/'^) is divided by the number of participant 
pairs Npartf^ in order correct for the small difference in the system size between Au + Au and 
Pb + Ph collisions. The data include fixed target measurements of Au + Au collisions at the 
AGS [SI at Brookhaven National Laboratory and Ph + Ph collisions at the SPS at CERN 
as well as colliding Au-beam collisions at four RHIC energies. First it is interesting to note 
that the measurements at ^smn = 200 GeV (the highest energy) falls in the lower part of the 
pre-2000 predictions (grey bars). This discrepancy between many of the predictions led to a 
strong revision of several model assumptions and discarded models which led to a severe over 
prediction of this simple observable. 

In addition, one observes that the energy dependence rather closely approximates a 
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Figure 3. The energy dependence of the mid-rapidity |t/| < 1 pseudo-rapidity density for 
central heavy-ion colhsions. Data obtained at RHIC (soUd circles) [HIS], SPS (open squares) [0], 
and AGS (open diamonds) |S] show a roughly logarithmic increase with collisions energy, y^sjviv 
(solid line) , but are also consistent with the trend of the Gluon Saturation Model (dashed curve) 
at y/WwN > 20 GeV. Corresponding data from pp ^W^ and pp ^2 collisions are shown as 
open symbols 



logarithmic increase denoted by the solid curve such that 

0.771n^/i^- 0.31. (2) 



2 C^A^c/i||r,|<l 



(Npart) dr] 

None of the present models reproduce this exact energy dependence. However, the gluon 
saturation model of Kharzeev, Levin, and Nardi |13[ I14j agrees rather well with the data above 
y/sNN > 20 GeV as indicated by the dashed curve in Fig. El This model leads to the expression 
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Here, Qs is the gluon saturation momentum, Aqcd is the QCD scale parameter, c is a 
normalization constant that is fixed at a single energy of y^=130 GeV, and A provides a 
scaling with energy, which for HERA data assumes a value of A ~ 0.25 — 0.3 The dashed 
curve shown in Fig. |21 corresponds to the following values; c = 1.02, A = 0.205 (obtained from 
the fit to the data), = 2.02(s/so)-^/^|lll, m = O.TSA/QlTc^^^^eV/c?, andpt = Qs- Although 
some of these values differ slightly from those listed in Ref . ) they are all considered to be in 
a reasonable range. 

In addition, one observes from Fig.lSJthat the particle production in heavy-ion collisions (solid 
symbols) is substantially higher that what is seen in pp ^01 orpp |lHll2j collisions (open symbols) 
when normalized to colliding nucleon-nucleon pair, Npart f^- One possible interpretation of 
this difference is that a large fraction (~50%) of the available energy is carried off by leading 
hadrons in pp and pp collisions, whereas this effect is absent in heavy-ion collisions because these 
hadrons suffer subsequent collisions leading to particle production in the heavy-ion collision 
environment jl5j 




Figure 4. The mid-rapidity density, dNch/drj scaled to the number of participant pairs, 
{Npartf^) is shown as a function of centrality, (Npart) for 200 GeV (left panel) and 19.6 GeV 
(right panel) j^. The solid squares (diamonds) represent the pp (pp) values at 200 (19.6) GeV, 
respectively jlll llOj . The solid curves show the scaling with the number of collisions obtained 
by the Glauber model using nucleon-nucleon cross sections of 42 mb and 33 mb for 200 and 19.6 
GeV, respectively 



3. Centrality dependence 

In heavy-ion collisions it is possible to obtain additional information on the particle production 
by varying the impact parameter between the two colliding ions, which also varies the size of 
the collision volume. In all four RHIC experiments, the centrality of a collision is measured on 
the basis of the total charged particle multiplicity or energy with in a certain pseudo-rapidity 
region. The methods for obtaining the centrality in the PHOBOS experiment is described in 
detail in a separate contribution to these proceedings [TBI and will not be discussed further here. 

As examples of the centrality dependence. Fig. shows the mid-rapidity density, 
dN / dr}\\^\^i/ {Npart 1'^) (solid points) as a function the number of participants, Npart for Au + Au 
collisions at 19.6 and 200 GeV [Hj. One observes a steady increase in the particle production 
with centrality in both cases, whereas scaling with the number of participants would lead to 
a flat dependence. This is shown by the dashed lines which correspond to the level seen in 
nucleon-nucleon collisions |^ (200 GeV pp), ^U] (19.6 GeV extrapolated pp). On the other 
hand, if the particle production would scale with the number of inelastic collisions between 
nucleons (obtained from Glauber model simulations), one would expect a very strong centrality 
dependence as shown by the solid curves in Fig. ^ Clearly, this strong dependence is not seen 
in the data; they are much closer to the Npart scaling limit. 

4. dN/drj Shapes 

The energy evolution of the shape of dN/di] distributions in central Au + An-collisions is 
illustrated in Fig. [3 Several observations can be made: 1) The mid-rapidity plateau increases 
in both width and height with collision energy and 2) the extend of the fall-off regions outside 
of this plateau also increases with energy, but the slope of the fall-off is essentially independent 
of energy, a consequence of limiting fragmentation scaling At first glance it may appear 
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Figure 5. The shapes of the pseudorapidity density, dNch/dr], is shown as a function of r] for 
0-6% central Au+Au colhsions at 19.6 (open squares), 130 (open circles), and 19.6 GeV (solid 
diamonds) .1,. The solid curves represent best fits to the data using triple Gaussians. 

strange that an order of magnitude increase in the collision energy results in only about a factor 
of two increase in the mid-rapidity value of dNch/dr]. Where does the energy go? Here it is 

worth keeping in mind that the energy of a particle is E = \J uiq + cosh^ rj, where mo and 
Pt is its rest mass and transverse momentum, respectively. The energy per particle therefore 
increases sharply with rj provided that {pt) does not fall off precipitously away from ry=0 . Thus, 
the energy of a pion with pt=0.5 GeV emitted at the half maximum point of ry ~ 2 at 19.6 
GeV is E=1.9 GeV whereas it requires an energy of E=11.2 GeV to emit the same particle 
at rj ~3.8, the half maximum at 200 GeV. The increased width of the dNch/dr] distribution 
therefore consumes the largest fraction of the additional energy. 

The evolution with both energy and centrality is shown in Fig. IS] for Au + Au collisions 
measured by PHOBOS (1^. Similar data at 130 and 200 GeV have been obtained by the 
BRAHMS collaboration [T71 UHl 

In Fig. 13 central dNch/dr] distributions at 19.6, 130 and 200 GeV are compared with model 
calculations. The dashed curves shown in the left panels are obtained from the gluon saturation 
model of Kharzeev et al. jl3| I14j using the analytical expression 



and the parameter values listed in Sect. 2. Note that a rigorous transformation from y to r] 
requires an integration over the actual pt distribution for each particle species, i.e. pions, kaons, 
and nucleons, which is not carried out here. We find that this analytical expression gives a 
reasonable account of the dNch/dr] distributions over most of the r] range, but that it deviates 
at the largest pseudorapidities, where the approximations on which it is based, are not fulfilled. 
In fact, a more accurate calculation by Kharzeev et ai. jl9j indicate that also the 19.6 GeV data 
can be reproduced to a satisfactory degree. 
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Figure 6. dNch/dr] vs r] (solid and open points) for five centrality bins representing 45% of 
the total cross section for ^sjvAr=19.6 GeV Au + Au collisions and for six centrality bins for 
y/sNN=130 and 200 GeV correcponding to 55% of the cross section p. The shaded bands 
represent estimated systematic errors. 



The 0-6% most central data are also compared to the HIJING 20 AMPT |2J models as 
shown in the right hand panels of Fig. EJ It is evident that the HIJING model (solid cm'ves) 
tend to underpredict the width of the distributions, whereas the AMPT model calculation at 
130 GeV appears to alleviate this discrepancy. Since the AMPT model uses the HIJING event 
generator but includes the final state interactions, it appears that this effect is important in the 
calculations. 

4-1. Boost invariance? 

Initially it was thought that a boost invariant region around mid-rapidity would develop at 
sufficiently high collision energy j7j. At first sight one might interpret the flat-top shape 
of dNch/drj distributions shown in Fig. \S\ as evidence for boost invariance, but this plateau 
is only a consequence of the transformation from rapidity to pseudurapidity space in which 
the Jacobian gives rise to a reduction of dNch/dr] near mid rapidity. This effect is clearly 
illustrated in Fig.|HI Here the dNch/dr] distribution for 0-6% central 200 GeV Au + Au collisions 
measured by PHOBOS T (open circles) is compared with the dN/dy distributions of pions 
and kaons measured by BRAHMS for 0-5% centrality for the same collision system. Whereas 
the dNch/dr] distribution exhibit a flat-top shape, a nearly Gaussian shape is observed for the 
dN/dy distribution, the width of which is in rather good agreement with the predictions of 
the hydrodynamical expansion model proposed by Landau |22j . In this model it is assumed 
that full stopping is achieved between the two colliding ions followed by isentropic expansion 
of a thermally equilibrated system. This leads to an approximately Gaussian shape of the 
dN/dy distribution [22|. In a simplified version of this model |2l] the standard deviation of 
the distribution is given by cj^ = \n{y/TNN /2mp), where rup is the proton mass. In the case of 
200 GeV Au+Au collisions, a = 2.16, the value used for calculating the solid curve in Fig.|HJ The 
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Figure 7. Pseudorapidity distributions of charged particles emitted in central (0-3% bin) 
Au+Au collisions at ^/sNN= 19.6, 130, and 200 GeV TH (symbols) are compared with saturation 
model calculations, Eqs. 7-9 (dashed curves). The solid curves represent multi-Gaussian fits to 
the data, see text. Pseudorapidity distributions of charged particles emitted in central (0-3% bin) 
Au+Au collisions at y/SNN= 19.6, 130, and 200 GeV (symbols) are compared with saturation 
model calculations, Eqs. 7-9 (dashed curves). The solid curves represent multi-Gaussian fits 
to the data, see text, are compared to calculations with the HIJING |20j (solid curves) and 
AMPT _21 (dashed curve) models. 

measured dN/dy distribution of pions is seen to deviate from this curve only slightly in the mid- 
rapidity region which may indicate a slight tendency towards boost invariance. Overall, however, 
the data show that the simplified picture of a wide boost invariant region at mid-rapidity does 
not appear to be valid. 

4-2. d+Au collisions 

In an attempt to understand the effects of secondary collisions in the Au + Au system, collisions 
between deuterons and Au ions have been studied at RHIC. It was necessary to use deuterons 
instead of protons in order to better match the mass to charge ratio of the fully stripped ^u-ions 
as both beams are bent in the same magnetic field near the interaction regions of the collider. 
In Fig. inithe dNch/dij distributions are shown for both minimum-bias d + Au collisions (open 
diamonds) and different centrality bins given as percentage of the total inelastic cross section . 
The method of determining the collision centrality is discussed in detail in Ref. ||1S|. Here the 
positive pseudorapidity corresponds to the deuteron beam direction whereas the Au ions go in 
the negative rj direction. It is evident that the most abundant particle production is found at 
negative rj values for both the minimum-bias distribution as well as for the distributions of all 





Figure 8. Visual comparison of the dN/dy distributions for pions (squares) and kaons 
(diamonds) measured by BRAHMS (8j with the dNch/drj distribution for all charged particles 
measured by PHOBOS P 0-5% and 0-6% central Au-|-Au colhsions at 200 GeV, respectively. 
Note the different axes used for the two different types of data. The solid curve represents a 
Gaussian with (7=2.16 (see text). 



but the most peripheral collisions. It is worth noting, however, that the reduction (or dip) near 
?7=0 is most likely just an effect of the transformation to pseudorapidity space; it is expected 
that this reduction is absent in dN/dy distributions for individual particle species. 

The asymmetric character of the dNch/dij distributions is a natural consequence of 
longitudinal momentum conservation in the collision. Thus for a central collision, one may expect 
that the two deuteron participants (one neutron and one proton) interacts with more than ten 
nucleons in the Au ion. The charged particles created in such a collision must therefore reflect the 
net momentum of the participants, which corresponds to p\\ = {N^^^ — N^^j.^) x 200GeV/c. This 
large excess of momentum in the Au direction therefore naturally leads to a more abundant 
emission of particles in the negative pseudorapidity region. The longitudinal momentum 
asymmetry is, however, strongly reduced for peripheral collisions where in the 80-100% bin 
we may expect that iVj^"^ ~ ^part ™d, indeed, we observe an almost symmetric dNch/drj 
distribution in this case (solid triangles in Fig. |^ . 

4-3. Shapes of "elementary" collisions 

In Fig.E3we compare the dNch/drj shapes of central (0-6%) Au+Au (solid circles) [T] and pp[TT] 
(solid diamonds) collisions at ^SiVAr =200 GeV. The Au + Au data have been normalized by the 
number of participant pairs {Npart/'^)- We observe that the additional particle production seen 
in Au + Au collisions (open circles) is peaked at rj=Q. It is thus reasonable to associate this 
additional production with secondary collisions in the Au + Au system because collisions tend 
to build up the transverse momentum component as the system progresses toward a thermal 
equilibrium. An overall increase of about 40% is seen in the Au + Au system. 

Although FigCni appears to indicate that the additional particle production in concentrated 
in the mid-rapidity region, the systematic errors associated with both measurements do also 




Figure 9. The dNch/dr] distributions for d + Au collisions at ^/Tmn =200 GeV as shown for 
different centrality bins j^. The grey bands represent statistical and systematic errors. The 
minimum bias distribution is shown as open diamonds. 




Figure 10. Comparison of the dNch/dri/{part/'^) distribution for 0-6% central Au+Au collisions 
aty/smf=200 GeV (solid circles) and 200 GeV pp fill (solid diamonds) collisions. The 
difference between the two distributions is shown as open circles. 
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Figure 11. Limiting fragmentation scaling for 0-6% central Au + Au jlj and pp/pp IIU) 
collisions. 



allow for an interpretation in which the additional particle production occurs over the while 
rapidity range with a constant enhancement factor of approximately 1.3 |15j. 

It has also been noted that leptonic collisions give rise to dNch/dyT distributions (where ut 
is the particle rapidity along the thrust axis of the outgoing particle jets) that are similar to 
those of pp/pp and heavy-ion collisions ^3]. In this comparison it should, however, be kept in 
mind that some distortions to the dN^h/dyT distribution will occur if they are converted into 
pseudorapidity space as illustrated above. 

J^.^. Limiting fragmentation scaling 

In a previous section we discussed the observation that the width of the dNch/dij distribution 
increases with collision energy. This effect is a consequence of limiting fragmentation scaling, 
according to which dN/dy distributions are identical in the fragmentation region close to the 
target/beam rapidity region, i.e., dN/d{y — ybeam) is energy independent p6^. Because for large 
?7 values y rj + lii{pt/mt), where the second term gives rise to only a small and essentially 
constant shift, one may also expect this scaling behaviour to apply to distributions measured in 
pseudorapidity space. This behavior has been observed first in pp collisions |llj . In Fig. Illh . 
we illustrate that this scaling extends to substantially lower energies by including pp data from 
Thome et al.^Q\. 

With the high-quality dNch/dr] data obtained by PHOBOS for Au + Au collisions it has been 
shown 1 that limiting fragmentation scaling also holds very accurately for heavy-ion collisions 
as illustrated for central collisions in Fig. Illb .b. The charged particle production thus increases 
toward mid-rapidity with the same rate of about a = d'^N^h/dr]'^ = 195 independent of collision 
energy. 

5. Total multiplicity 

As observed on Figs. El and El the acceptance of the PHOBOS setup covers essentially the full 
distribution of charged particles emitted from Au + Au collisions, even at the highest energy 
of y^SiVAr=200 GeV. The unobserved fraction may be estimated in several ways ^^l; namely by 
integrating fitted Woods-Saxon curves or by making use of the limiting fragmentation scaling 
to estimate the tails of the distributions. Several such methods give almost identical results. 
Thus it is found that the extrapolated region accounts for less than 1% of the total for the most 
central bin increasing to < 5% for peripheral collisions at 130 and 200 GeV. For 19.6 GeV it is 
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Figure 12. Upper panels: The total charged particle multiplicity |15j (solid points) is shown as 
a function of centrality, {Npart) for three energies. Middle panels: Same as above but normalized 
to {Npart) f^- Lower panels: Centrality dependence of the full width at half maximum (FWHM) 
of the dNch/drj distributions (solid points). Solid curves represent HIJING calculations 201) 
whereas solid squares correspond to pp ^] collisions. Solid diamonds represent an extrapolation 
to 19.6 GeV oipp data UHl. 



believed that 100% of the charged particles originating from the collision region fall within the 
detector acceptance. In this case the total number of charged particles, Nch-, was estimated by 
integrating over the PHOBOS acceptance, i.e. |r/| < 5.4. It is however, believed that the tails 
seen at large pseudorapidities, especially in peripheral collisions, are associated with emission 
from the spectators. This method may therefore lead to a slight over estimate of Nch for the 
most peripheral collisions. 

The results are summarized in Fig. El where the total number of charged particles, Nch, 
is shown as a function of centrality expressed in terms of {Npart) in tbe three upper panels 
(solid points) |15j . The data are compared to HIJING calculations (solid curve) and to pp (solid 
square) and an extrapolation of pp (solid diamond) JHI collisions. By normalyzing to the 
average number of participants pairs {Npart 1"^) an interesting trend becomes evident; unlike the 
mid-rapidity density per participant pair dN / dri\\ri\<i/ {^part/"^) shown in Fig.|3J which increases 
substantially with centrality, we find that the total normalized multiplicity remains constant with 
centrality within experimental errors. This requires that the width of the distributions become 
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Figure 13. Nch/ {Npartf^) for Au + Au TF (solid points), d + Au j2Sl (open circles) are shown 
as a function of (Npart) and compared to valued for pp (TTI (solid square) and e+e" [27j (solid 
diamonds) collisions. The horizontal dashed line represents the average value for Au + Au 
collisions, which is ~40% higher than the pp level (solid line). 



narrower with centrality in order to keep the total area constant. This trend is indeed observed 
in the bottom three panels of Fig. ^Jwhich show that the full width at half maximum, FWHM, 
of the distributions [T] decrease with cenrality. It is interesting to note that the constant value 
of Nch/ {Npart 1'^) occurs at a substantially higher level than the pp data point, and that the 
HIJING model does not predict this constancy with centrality. 

It is also of interest to determine whether the total charged particle production in the d + Au 
collisions shows the enhancement over simple nucleon-nucleon collisions seen for heavy ion 
collisions. The total number of charged particles is in this case estimated by an extrapolation 
outside the measured region (shown in Fig. |^ based on the limiting fragmentation scaling using 
lower energy measurements of proton induced collisions j25j. The total normalized multiplicity, 
^ch/ {NpO'i^t f^) for d + Au is in Fig. [T^ shown as a function (Npart) (open points) and compared 
to the data for Au + Au (solid points), pp 11 (solid square) and e~^e~ j2Zl (solid diamond) 
collisions. We observe that the also the d + Au multiplicity is substantially lower than for the 
heavy-ion system as was seen in the comparison to pp collisions. This effect does not appear 
to be simply a consequence of the number of participants in the collisions because there in no 
indication of an increased multiplicity even for the most central d + Au collisions. We also note 
that the multiplicity for e~^e~ collisions is consistent with the heavy- ion data. As mentioned 
earlier, one possible explanation for this disparity between nucleon-induced and heavy-ion or 
lepton induced collisions is the fact that leading particles in nucleon-induced collisions carry 
away about half of the energy, which is then not available for particle production JHI- It is 
surprising, however, that this mechanism is still effective in d + Au collisions, which involve a 
substantial number of subsequent collisions. 

Earlier, we have noted that shapes of the dNch/dij distributions are to a reasonable 
approximation represented by a mid-rapidity plateau followed by a nearly linear fall-off in the 
fragmentation region. It is therefore tempting to approximate this shape by a trapezoidal 
distribution as illustrated in Fig. We have also seen that the mid-rapidity density 

dNch/d7]\^^^^i/ {Npart Z"^) increases logarithmically with ^/sJm (see Fig. such that the height 
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Figure 14. dNch/dij for 0-6% central 200 GeV Au + Au collisions T (open circles) is seen to 
be well approximated by a trapezoidal shape (solid lines) for which the parameters are given in 
the figure. 



of the midrapidity plateau increases logarithmically with energy (see Eqn.l^J, which corresponds 
to dNch\\Ti\<i/dr] = 132 In ^/snn — 53 for 0-6% central Au + Au collisions with Npart='iAA. In 



addition, the linear part of the limiting fragmentation region is well represented by 

2 dNch 



{Npart) dr) 



1.134(y6eam +0.15-7?), 



(5) 



which for 0-6% central Au + Au collisions corresponds to dNch/dr] = I95{ybeam + 0.15 — r/) as 
illustrated in Fig. El The area of the trapezoid, i.e. Nch, may then be computed as 



which corresponds to 
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for 0-6% central Au + Au collisions. 

As discussed in Sect. 2, the gluon saturation model 13 also gives a good representation of 

the mid-rapidity density ^jv^ f) '^^"^^'''^^ at energies above y^iAw=20 GeV. Using this estimate 
instead we obtain the dashed curve in Fig. El which also gives a good account of the total 
multiplicity. In conclusion we find that the simple trapezoidal representation of the dN^h/dr] 
distributions in heavy-ion collisions gives a good representation of the energy evolution of total 
charged particle production. 




Figure 15. Nch/ {Npart/'^) for central Au + Au and Pb + Ph collisions is seen to be well 
approximated by the expression given in Eqn. IHl (solid curve). The dashed curve corresponds to 
the gluon saturation model estimate of dNch/dij. 



6. Summary and conclusion 

In this work we have attempted to give a summary of the present status of charged 
particle multiplicity measurements in heavy-ion collisions with special emphasis on the recent 
measurements obtained at RHIC, but also studying how these results compare with lower energy 
measurements of Au + Au and Pb + Pb collisions studied at the AGS and SPS, respectively. 
Parallels and differences between the heavy-ion data and those obtained in nucleon-nucleon, 
nucleus-nucleon, and lepton-lepton collisions have been discussed. The first result from the new 
RHIC facility indicated a logarithmic trend in the energy evolution of the mid-rapidity density 
that was not predicted by most theoretical models. This logarithmic increase with collision 
energy has been verified in subsequent measurements and has helped to focus the attention 
on models that approximately reproduce this trend, one such model being based on the gluon 
saturation concept I14| I19j . The mid-rapidity density in Au + Au collisions shows only a 
rather weak centrality dependence when normalized to the number of participant pairs Npartf^', 
for the most central collision we observe at y^SAr/v"=200 GeV a 40% enhancement over that 
for nucleon-nucleon collisions at the same energy. The observed dependence is, however, much 
weaker than a scaling with the number of nucleon-nucleon collisions would suggest. 

The distributions of charged particles in pseudo-rapidity space exhibit a midrapidity plateau 
followed by an almost linear fall-off to higher values of \r]\ corresponding to the fragmentation 
region. A comparison with distributions in rapidity space of identified pions and kaons performed 
by the BRAHMS collaboration show that the mid-rapidity plateau seen in //-space should not 
be interpreted a evidence for a boost invariant region. 

A comparison of dNch/drj distributions at three different energies in Au + Au collisions 
demonstrate that the limiting fragmentation scaling observed earlier in nucleon-nucleon collisions 
hold rather rigorously also in heavy-ion collisions. Measurements of heavy-ion dNch/dr] 
distributions allow for a small extrapolation to obtain the total number of charged particles 
emitted from the collision region. It is found that for 200 GeV Au + Au this quantity, scaled by 
the number of participant pairs, Npartf^, is essentially constant as a function of centrality, but 
at a level of about 40% higher than for nucleon-nucleon and deuteron-gold collisions. Finally, 
we propose a simple expression, which is based on the observed approximate trapezoidal shape 



of the dNch/dr] distributions in Au + Au collisions, for which the height of the mid-rapidity 
plateau is given by the observed logarithmic dependence on the collision energy and the width 
is determined by the limited fragmentation scaling. This phenomenological expression provides 
an excellent account of the total multiplicity in heavy-ion collision over two orders of magnitude 
in collision energy. 
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